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I. INTRODUCTION 



The problem of test particle diffusion in stochastic magnetic fields was studied by many authors [|13[ and im- 
portant progress was obtained. However, the general solution was not yet found. Particle trajectories in a magnetized 
plasma are determined by three stochastic processes: the magnetic field, the collisional velocity along magnetic lines 
and the collisional velocity perpendicular to the magnetic lines. These components of the stochastic collisional ve- 
locities have very different effects. There are two important difficulties appearing in this triple stochastic process. 
One is related to the parallel collisional velocity which enters as a multiplicative noise in the equations of motion and 
the other with the Lagrangian non-linearity which is determined by the space dependence of the stochastic magnetic 
field. Each of these two problems has been recently studied but only considered separately. The complete model for 
particle transport in stochastic magnetic fields could not be analyzed until now. 

The latter difficulty can be eliminated if one restricts the study to stochastic magnetic fields with small amplitudes 
and/or large perpendicular correlation lengths for which the magnetic Kubo number (defined below) is small. If the 
perpendicular collisional velocity is neglected, this quasilinear problem has an exact solution that was obtained by 
several methods [Q. It shows that the parallel collisional motion determines a subdiffusive transport across the 
confining magnetic field with the running diffusion coefficient D{t) decaying to zero as D{t) ~ t^^/^. It was shown 
p3[ that this subdiffusive transport is due to collision induced trajectory trapping along the magnetic lines. The 
parallel collisional velocity forces the particles to return in the already visited positions along the magnetic lines 
and consequently generates long time Lagrangian correlation of the stochastic magnetic field. If the perpendicular 
collisional velocity is taken into account, the transport is diffusive and the diffusion coefficient was evaluated semi- 
qualitatively by several methods Q. 

On the other hand, the Lagrangian non-linearity determined by the space-dependence of the stochastic magnetic 
field leads, at large magnetic Kubo numbers, to magnetic line trapping and generation of stochastic magnetic islands. 
This process is mathematically identical with the trajectory trapping in the E x B drift motion in an electrostatic 
turbulence. The latter was recently studied by means of a new statistical approach, the decorrelation trajectory 
method [0. 

The aim of this paper is to study the general problem of collisional particle diffusion in stochastic magnetic fields in 
the guiding center approximation. More specifically, we determine the effect of self-consistent generation of stochastic 
magnetic islands on the effective transport. The running diffusion coefficient is determined for arbitrary parameters 
of the stochastic magnetic field and of particle collisions. The decorrelation trajectory method is used for studying 
this rather complicated triple stochastic process. 

The paper is organized as follows. The physical model is described in Section 2. We derive in Section 3 the 
Lagrangian velocity correlation and the running diffusion coefficient for arbitrary values of the four specific parameters 
and for given Eulerian correlation of the potential. The physical significance of this general result is then analyzed: 
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the subdiffusive transport in Section 4, the effect of coUisional cross-field diffusion in Section 5 and the effect of a time 
variation of the stochastic magnetic field in Section 6. The conclusions are summarized in Section 7. 



II. THE SYSTEM OF EQUATIONS 

The particle guiding center motion is studied in a magnetic field with a stochastic component. The magnetic field is 
taken to be a sum of a large constant field Bo=-Boez and a small fluctuating field perpendicular to Bq, and depending 
on the perpendicular coordinates x = (a;, y) and on the parallel coordinate z 

B = Bo(ez + b(x,0,i)) (1) 

(Here the perpendicular and the parallel directions are defined in relation to the direction of Bq). This is the usual 
slab model of the confining c;onfiguration in a tokamak plasma. Since the reduced magnetic field is divergence-free, 
V • b = 0, its two components can be determined from a scalar function (/)(x, z) as 

b(x, = V X (j!)(x, 2;,i)e2. (2) 

The system of equations for guiding center motion is: 

dx 

— =b(x,z,f)ry||(i)+,?i(t), (3) 
dz , , 

^=.||(*). (4) 

The three stochastic functions b(x, z, t), ri±{t) and ryu (t) are statistically independent: all cross correlations are zero. 
All these stochastic fimctions arc assumed to be Gaussian, stationary and homogeneous, with zero averages. The 
autocorrelation function of the stochastic potential (/)(x, z,t) is modeled by: 

^(x, z, t) = (</>(0, 0, 0) 0(x, z, t)) = (3^X1 exp (^"^ " ^^"^TTT^ I exp ( ) (5) 

where f3 is the mean sqiiarc value of the reduced magnetic field b, Ay is the correlation length of the potential (p along 
the main magnetic field Bq, Ax is the correlation length in the plane perpendicular to Bq and Tc is the correlation time 
of (p. The autocorrelation tensor of the reduced magnetic field components Bij = {bi{O,0,0)bj{x,z,t)) , i,j = x,y, is 
determined from A{x, z) as 

g2 Q2 q2 

The coUisional velocities are modeled by colored noises with the correlations 

{v\\{0)V\\{t)),^X\\'^R{'^t) (7) 




(^vi{0)rj{{t))^ = 6ijx±uR{ut) (8) 

where i' is the collision frequency, X|| = ^mfp'^/'^ parallel coUisional diffusivity, Xmfp is the parallel mean free 

path, x_L = p|i//2 is the perpendicular coUisional diffusivity and pi, is the Larmor radius relative to the reference 
field. R{i't) is a time decreasing function that is chosen as 

R{ut)=exp{-i/\t\) (9) 

for the explicit calculations presented in this paper. 

We introduce dimcnsionless quantities with the following units: A_l for the perpendicular displacements, A|| for the 
displacements along the reference magnetic field and for the time. The perpendicular velocity v = br/|| is reduced 
with V = /3y/xj|J7, the parallel velocity r^y (t) with ^/XjfJ^ and the perpendicular coUisional velocity f]± {t) with ^/xii'. 
The equations of motion in these dimensionless variables (denoted by the same symbols as the physical ones) are: 
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^ = Mb(x,z,%|(i)+xi/VW (10) 

§-xJ^'^iiW- (11) 

Four dimensionless parameters appear naturally in this problem: the dimensionless perpendicular and respectively 
parallel diffusivities 

a dimensionless parameter that contains the effect of the stochastic magnetic field 
and the dimensionless decorrelation time: 

fc = T^V. (14) 

We note that the parameter which describes the evolution of the magnetic lines, the magnetic Kubo number = 

1 /2 

/3A[|/Aj^, appears here as a factor in M, which can be written as M — KmX\\ ■ 

The aim of this calculation is to determine the Lagrangian correlation of the effective perpendicular velocity 

v(x, z,i) = b(x, z,i)?7||(t) (15) 

which leads to the perpendicular effective diffusion coefficient. 

III. SOLUTION BY THE DECORRELATION TRAJECTORY METHOD 

We use the decorrelation trajectory method following the recent calculations for the influence of particle collisions 
on the diffusion in electrostatic turbulence [^. The difference and the supplementary difficulty of the magnetic 
problem comes from the structure ( p^ of the velocity v which is the product of two stochastic processes. They are 
statistically independent but in the Lagrangian frame they are correlated through the trajectories due to the space 
dependence of the magnetic field fluctuations. The later makes this problem strongly nonlinear. The trajectories also 
depend on the coUisional velocity r]± and thus the velocity v is a triple stochastic process in the Lagrangian frame. 

We determine the collisional contributions to the perpendicular displacement: 

m=xT l\dr)dT (16) 



and make the change of variable x'(t) = x(t) — ^{t) in Eq.(lO), which introduces the collisional displacements in the 
argument of the magnetic field fluctuations: 

— = M^[^[t)+m.zA- (17) 

Here v [x'(i) + ^(i), z,t] — h [x'(i) + £,[t)^ z, t] Tyy (i) is a triple stochastic process. 

We calculate first the Eulerian correlation (EC) of b [x + ^(t), z, t] which is defined as an average over the magnetic 

field fluctuations and over the perpendicular collisional velocity. We calculate the EC of the potential 0(x, z,<) = 
[x + C(i)i A ^'^d then derive the EC of the magnetic field components. 

E = ((0 [Xi + ah), Zi, <i] Ct, [X2 + C(t2), Z2,t2]))c = 

= {A [Xi + e(il) - X2 - i{t2), Zi - Z2, h - t2]), (18) 

This average over the perpendicular collisional velocity can be calculated using the 2-point Gaussian probability 
density: 
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E= / deide2^[xi- x2+a- 6,^1- Z2,ti-i2] ^'2(^1,^1; 6,^2) (19) 



where -P2(Cij C2, ^2) is the probabihty density for having £,{ti) — £_i and ^(^2) = ^2- It is dctermmed as the average 
over colhsions of the corresponding product of 5-functions: 

This probabihty can be calculated using the Fourier representation of the ^-functions and the cumulant expansion of 
the resulting exponential. Since the coUisional displacements are Gaussian, only the first two cumulants appear. One 
obtains the two-point probability density for the perpendicular coUisional displacements as 

P2 - / / dqidq2 exp ( *qi • 6 + *q2 • 6 - ^^^'^^^ " ^^^'^^^ " (qi • ^(ii) q2 • ^(^2)) J (20) 



which, introduced in Eq.(19), yields 

i?(x, z, r) = y de ^ (x + e, z, t) Px(C, r) (21) 

after using the dependence of the EC on the difference ^ = £.1 — ^2 and performing the integrals over qi, q2 and ^1. 
Here P±{£,,t) is the one-point probability density for the perpendicular coUisional displacements: 



2^(C2(r))^ 2(e^(r)) 
The MSD for the coUisional perpendicular displacements is: 

(K(t2) - mf)^ = (eir))^ = X^*(r) (23) 

where r = 1^2 — ii| and ^'(t), the reduced mean square coUisional displacement, is: 

^(t) = 2 /" (r - t)R{t)dt = 2 [t + exp(-T) - 1] 
Jo 



(24) 



Thus the average effect of the perpendicular coUisional velocity ri^{t) consists of the modification of the EC of 
the magnetic potential (j). The EC j4(x, z,r) is transformed into E{x.,z,t) [Eq.(|2l])] gaining a supplementary time- 
dependence in addition to the one determined by the finite correlation time of the stochastic magnetic field. As observed 
in ||l^ , E is the solution of a diffusive equation and the effect of collisions consists in progressively smoothing out 
the EC of the magnetic potential and in eliminating asymptotically the x de pen dence of i?(x, z, r). Since the integral 
over X of -E is constant, the time dependence introduced by collisions in Eq.(21) does not destroy the correlation but 
only spreads it out. 

We note that the average over the coUisional parallel velocity was not performed at this stage: z is in Eq. (|2^) an 
Eulerian coordinate. 

The problem of coUisional particle motion in magnetic turbulence (|To[), (|ll|) is now formally reduced to a doubly 
stochastic process: 

dx ~ 

— = A/b(x,2,i)77,i(i) (25) 



i-xi^llW (26) 

where h{yi,z,t) is the stochastic magnetic field generated by the potential 0(x, z,t). The effect of the perpendicular 
coUisional velocity is an additional time-dependence introduced in 0(x, z,t) and the transformation of its EC from 
Eq.(|l8|) into Eq.(|2^). The Eulerian correlation of the components of b(x, z,i) are determined from the EC of the 
potential ( ^ ) by equations similar to (H). 

The Langevin equation ( p5| ) can be written as dx/dt = v(x, t) and thus it is similar to the two-dimensional 
divergence- free problem studied in . The velocity 
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v(x, z, t) = b(x, z, t)77|| (t) 



(27) 



has a much more comphcated structure being determined by two muhiphed stochastic processes. However, the method 
developed in can be used here: we wiU follow the same calculation steps as in 

First, we define a set of subensembles S of the realizations of the stochastic functions that have given values of the 
potential 0, of the magnetic field b and of the parallel velocity r/|| in the point x = 0, z = at time t — : 



(0,0,0) = 0", b(O,0,0)=b", miO)^V°- 



(28) 



The correlation of the Lagrangian velocity (g^) can be represented by a sum over the subensembles of the correlations 
appearing in each subensemble 



L{t) ^ J d0"db"d?7°P(b°, 0^ 77") (v(0, 0, 0)v [x{t), z{t),t] 



(29) 



where P (b", 0°, 77") P (bf) P (6^) P (0°) P (77°) with P {X) = exp {-X^/2) /\/2^ is the probabihty of having 
b°,(/)°,77" at X = 0, z = and t = 0. This probability is a product of individual distributions because the stochastic 
variables are not correlated in x = 0, z = 0, t = 0. The point x = 0, z = is taken as the initial condition for the 
trajectories determined from Eqs. (psf), (p6|). Since the initial velocity in the subensemble S is v(0, 0,0) = b°77° for 
all trajectories, the subensemble average in Eq.(|29|) is (v(0, 0, 0)v(x(i), z(t), t))^ = b°77° (v [x(t), z(t), and thus 
the Lagrangian correlation L(t) is determined by the average Lagrangian velocities in all subensembles. In order to 
evaluate these quantities, we need to calculate the average Eulerian velocity in the subensemble S, 



V5(x, t) = (v [x, z{t),t])s - (b [x, z{t),t] 7711 (t)) . 



(30) 



where is the average over the two stochastic processes restricted to the realizations in S and z(t) is the stochastic 
parallel displacement obtained from Eq.(p^) 



dTm{T). 



(31) 



More precisely, Vf is determined by the following conditional average 



b,[x,z{t),t]r]\\{t)5 bO-b(0,0,0) 6 - cj) (0,0,0) S [jf - ^ 



P(bO,(/)0,770) 

Introducing a function S [z — z(t)] and using the statistical independence of b and 771 one can write 



(32) 



vf ^ dz (6i(x, z, t) 6 \h° - b(0, 0, 0) 



(0,0,0) 



1 



P(bo,0") 



(7711 (t) S [z - z{t)] 5 [77" - 7711 (0)] )^ (33) 
The first average over the stochastic magnetic field represents the subensemble average of 6i(x, z,t) in S and is given 

by 

Bf{x,z,t)={h{x,z,t))g=(^, coS(x,z,t) (34) 



dx2 ' dx 



where $'^(x, z,t) = ((/)(x, z,t)j , the average potential in the subensemble S, is calculated as in reference and is 
the following function of the parameters 0", b° of the subensemble and of the EC of the potential (pi|): 



where Ei^{x,z,t) = (6^(0, 0, 0) (/)(x, z, t) ) ~eij-^E{x, z,t) 



$^(x, z, t) = (t>°E{ji, z, t) + h\E,^{x, z, t) (35) 
a 
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The second average in Eq. (^3|) over the coUisional parallel velocity can be written using the Fourier representation 
of the (5-functions as 



M„ ^{,711 (t)<5[z-z(t)]<5 [r?0- ,7,1(0]), 



1 



P(?7") 



dkdq exp [—ikz — iqnf^) (?7j| [t) exp \ikz{t) + iqrj\\ (0)] 



(36) 



The average in this equation can be calculated as the derivative with respect to a of the following average, evaluated 
in a = 

(exp (a77|| {t) + ikz{t) + iqri\\ (0)))^ = 



where 



exp 



o?" k'^ 

-y - Y W)c - y + (m {t)z{t))^ + iaqRit) - kq (77,1 (O)z(i)), 



(^||(0)^(t)),-xi/' / dTR{r)^x\Cv{t) 



,1/2. 



(37) 



(38) 



(r7|| {t)z{t))^ = x\'^ I drR{t - r) = Xy^'^^W 



,1/2. 



(39) 



{z\t))^= dndT2R{\T^-T^\)=X\\^(t). (40) 

Jo 

For the correlation R in Eq.(|^), the reduced parallel running diffusion coefficient is 

V{t) = 1 - exp(-t) 

and the reduced mean square parallel displacement is ^'(t) defined in Eq.(|2^), the same as for the perpendicular 
coUisional displacement. 

Straightforward calculations lead to the following equation for the parallel average Mii |3^: 



M||(z,<) = 



,^Rit)-x];Mt)l^ 



pS{z,t) 



(41) 



where P^{z,t) is the probability of having a parallel displacement z at time t taken for the trajectories in the 
subensemble S. This was obtained as a Gaussian distribution with an average displacement {z{t)) g and a modified 

/ 2\^/2 

dispersion s{t) — {^{z{t) — {z{t))g) y : 



P'iz,t) 



1 



■ exp 



The parallel average displacement is the integral of the parallel average velocity in S 

{m{t))g = ,fR{t) 

and is obtained as 

The mean square parallel displacement of the trajectories in S is 



(42) 

(43) 
(44) 
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\t) = {z'{t)) - Xl\'D\t) ^ X|| (*(r) - V^{t)) . 



(45) 



Thus the dispersion of the paraUel component of the trajectories in a subensemble S is always smaller than the 
dispersion of the whole set of trajectories <^z^(t)). It grows slowly (as t^) at small t and at i ^ 1 it reaches (z^(i)) . 
The parallel running diffusion cocfhcient in the subensemble S is 'D^^{t) — X||2-'(i) (1 — It behaves at small time 

as and at t 1 it is equal to the global diffusion coefficient of the whole set of trajectories. 
The average velocity ( p2| ) in the subensemble S is thus obtained using Eqs. (|33|)-(^l|) as 

(46) 



l/,^(x,t) = J dzBf{^,z,t)Amz,t). 

The next step in the decorrelation trajectory method is to find a deterministic trajectory X'^(t) in each subensemble 
S as the solution of the equation 



dt 



MV^(X'^,t) 



(47) 



with X'5(0) = 0. Using Eqs.(g|) and (g) one can show that this is a Hamiltonian system of equations which can be 
written as: 



dt 

dY^ 
dt 



-M 
M 



9g^(x^,y^,t) 



(48) 



with the Hamiltonian 



H^{X^,t) = J dz$^(X^,z,t)M||(z,i). 



(49) 



This Hamiltonian represents the average potential in the subensemble S. Its explicit expression calculated for the 
correlations (||) and (^) is: 



H^{X'^,y'^,t) = br^°{p-n^Y'^)n^exp[--n^[{X'^) + (F^) ) /|| 



(50) 



where 



/||(t) =nJ/'(i?-Xp||I?2(i„i?))exp -i (,;")' 



(51) 



n^it) = [l + x^^{t)] \ n||(i) = [l + s2(t)]- 



(52) 



Since the stochastic magnetic field considered here is isotropic, the Hamiltonian could be simplified by taking the x 
axis along b*'. The parameters of the subensemble S are in Eg. (50) b = \h^ \ , p = cfP/b and 77°. The equations for the 
decorrelation trajectories (^ ) obtain from the Hamiltonian (js^) are 



dX^ 

dt 
dYS 

dt 



X 



= -MbrPniS\\X^ (p - 7ij_r^) exp ( --n^ \^{X^) + (F ' 



(53) 
(54) 



The average Lagrangian velocity is estimated as in p6| by the average Eulerian velocity along the decorrelation 
trajectory 



{^[^{t),t]), = Y' [X^(t),t] 
where X^{t) is the solution of Eqs.(|3|)-(||). 



(55) 
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We finally obtain using Eq.(55) and ( ]29| ) the correlation of the perpendicular Lagrangian velocity v [x(t), z{t),t] ri\\ (t) 
for arbitrary values of the four dimensionless parameters (12)-(^^ and for given Eulerian correlations of the three 
stochastic processes that combine in the equations of motion @-^): 

L{t;M,X\\,X±,Tc) = {i^X±fM^^x 



27r 



dp I db exp ( -y (p^ 



drprp exp 



„02 



(56) 



The total perpendicular running diffusion coefficient is the sum of two terms: a direct contribution of the collisional 
velocity r]±_ obtained from Eq.(p3|) and the contribution of the velocity (p^: 

D{t- M, x\i,x±,rc)^ XiT^{t) + {y\\)M D,nt {t; M, X|| , X± , Tc). (57) 

The latter is the time-integral of the Lagrangian correlation (|5^) and can be written as: 

1 . / b^ 



= ^ y dp J d6 6^exp |^-y(p2 



1 



difrp exp 



(58) 



where X^{t) is the component along x axis of the solution of Eq.(^). It depends on the parameters M, x\\i X± a-nd 
as well as on the shape of the Eulerian correlations. This contribution ( |58|) results from the nonlinear interaction 



of the three stochastic processes. These results (|56|)-(57) are written as dimensional quantities. 

A computer code that calculates the running diffusion coefficient starting from the analytical expression ( |58|) was 
developed. It determines the decorrelation trajectories (48) for a large enough number of subensembles and performs 
the integrals in Eq.(58). The code was tested and the parameters in the numerical calculation were established using 
the analytical results concerning the subdiffusive transport. Namely, as shown in the next section, the asymptotic 
expression for the decorrelation trajectories and for the diffusion coefficient can be determined for arbitrary M and 
X\\ if X_L = and Tc = oo. This provides a very good test for the code and permits the optimization of the choice of 
the parameters. 

The analyses of the collisional particle transport in stochastic magnetic fields obtained by means of the decorrelation 
trajectory method results (pfl)-(58) is the subject of the next three sections. 



IV. SUBDIFFUSIVE TRANSPORT 



We first consider a static stochastic magnetic field (tc — > oo) and the zero Larmor radius limit corresponding to 
negligible cross field collisional diffusion, xi- = 0- It is interesting to study separately this particular case because it 
leads to a subdiffusive transport determined, as shown below, by two kinds of trapping processes. Moreover, the time 
dependence of the diffusion coefficient obtained for these particular conditions allows the understanding of the scaling 
lows of the diffusion coefficient determined by the presence of a decorrelation mechanism. 

For the limit Aj_ — > oo, an exact analytical solution was determined fl^ . It was shown that particle perpendicular 
transport is subdiffusive with the running diffusion coefficient going asymptotically to zero as t^^/"^. This particular 
case is used here as a test for the decorrelation trajectory method. We show that the exact solution is found. Then 
the non-linear problem corresponding to finite A_l is studied. We show that the generation of magnetic islands by 
magnetic line trapping does not change the asymptotic behavior of the diffusion coefficient: a similar subdiffusive 
regime is obtained with D{t) ~ f^l"^ . The nonlinear process of island generation has a strong effect but it is localized 
in time: it determines a transient decrease of D(t). This effect is very important because it leads, as will be shown 
in the next sections, to complex anomalous regimes when 7^ or when is finite. 

In the limit Aj_ — > 00 the Lagrangian non-linearity determined by the x-dependence of the stochastic magnetic field 
disappears and the problem simplifies considerably. The equations for the decorrelation trajectories (^^ reduce to 

f = -^-»/„(0. f = (59) 

where dimensional quantities were used. Thus the average Lagrangian velocity in S needed for determining the 
Lagrangian velocity correlation according to ( p6| ) is Vf^(t) = — 6r/°/||(i). The integrals over p, b and can easily be 
performed in Eq.(pq) and one obtains 
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Loit; 0, X|| , 0, oo) = ; 

[l + X||*W 

which after algebraic transformations becomes 



3/2 



(60) 



Lo(i;0,X||,0,oo) = V'- 



1 + X||*W 



1/2 



i + X||*W 



(61) 



This is precisely identical with the exact analytical solution determined in |14| by means of a different method. The 
perpendicular running diffusion coefficient can be obtained by time- integration of Eq. ( |6l| ) as 



1 + X||*(^) 



1/2- 



(62) 



1 /2 

This exact solution obtained for Aj^ — > cx) is also valid for finite Aj^ if M = /3A||/Aj^X|| <C 1. Actually this is the 
condition for neglecting the perpendicular displacements and the x-dependence of the magnetic field fluctuations. 
Consequently, Eqs.(|6l|), (|6^ ) have physical relevance for tokamak plasmas, although Aj_ is of the order of 1 cm and 
it is smaller than A|| by at least a factor 10'^. Due to the small values of (3 which are usually of the order 10"^ the 
parameter M can be small. 

The absolute value of Lo(^; X|| , 0, oo) and DQ{t; M,x\\,0,oo) are plotted in Figures 1 and 2. One can see that 
the Lagrangian correlation has a long negative tail at large t; its contribution exactly compensates the positive part 
appearing at small time such that its time-integral is zero. More precisely, Dq ~ for long time. The zero of the 

Lagrangian correlation (and the maximum of Do) appears at the average returning time r^. It is determined from the 

— 1/2 

equation Lo(rr ; M, x\\,0,oo) = and it is a decreasing function of X\\ scaling approximately as X|| . It is remarkable 
to note that in the limiting case of absence of collisions (i^ = 0), Eq.(|62[) yields a finite diffusion coefficient. In this 
case, x\\^ — (where Vt is the thermal velocity) and a small time expansion can be done in Eq.(|62[) obtaining 

the result of Jokipii and Parker Jp], Djp = /3^A||Vt/V2- This is also well known as the Rochester and Rosenbluth 
collisionless diffusion coefficient [Qf , in the form Drr — DmVr, where Dm is the diffusion coefficient of the magnetic 
lines (see also |l^, ||l^]). Thus, the collisions determine a very strong change of the perpendicular transport, which 
is diffusive in the absence of collisions and becomes subdiffusive due to the parallel coUisional motion. A physical 
interpretation of this subdiffusive behavior is presented in in terms of a parallel trapping process determined by 
the collisions which force the particles to return in the already visited points along the magnetic lines. Consequently 
the Lagrangian velocities remain correlated. Since the parallel velocity changes its direction due to collisions, this 
long-time correlation is negative and thus determines the decay of the running diffusion coefficient Dq (t) . 

A similar subdiffusive transport appears in the non-linear case too, provided that x± = 0- f^is case, n±{t) — 1 
and the Hamiltonian ( pO| ) depends on time only through the factor /||(<). It can be written as: 



H{X,Y,t) ^ f^{X,Y) /||(t), 
and consequently one can make a change of variable from t to r(i) defined by 

and the equations for the decorrelation trajectories become: 



dX 
dT 



-M 



df±{X,Y) 
dV 



dY 
dT 



M 



df±iX,Y) 
dX 



(63) 



(64) 



(65) 



The function T{t) has a maximum and then decays to zero. The solution of the time-independent Hamiltonian 
equations ( |65| ) is a periodic function of r with X'^(t) lying on the closed paths determined by f^(X,Y) = /^(0,0). 
The size of the paths depends only on p : it is infinite (straight line) &t p = Q and decays to zero as p increases. 
The period is proportional to [Mhrf )^^ . The decorrelation trajectories are thus obtained as X'^(T(t)) where X'^(t) 
is the solution of (|65|). This show that the trajectories wind around the closed paths (for an incomplete turn or 
for many turns, depending on M and on the parameters p, b and 7]^); at the time corresponding to the maximum 
of T{t) they all stop and go back along the same path. Since r(<) when t oo, the asymptotic value of the 
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abs[L] 




FIG. 1. The Lagrangian velocity correlation for the subdifFusive transport ( X\\ = Q, Tc = oo). Lo(i) corresponds to M <^ 1 
and is given byEq.(|ll) and is the nonlinear result obtained in the presence of magnetic line trapping at large Km {M = 10, 
X|| ~ 0.1). The dashed parts of the two curves represent negative values of the Lagrangian correlations. 



decorrelation trajectories is X'^(T(i)) — > X'^(O) — 0. All decorrelation trajectories eventually stop at the origin. The 
equation for the diffusion coefficient ( p7| ) thus gives D{t) 0. Using Eqs.(|6^) and (|T]) the function r(t) is shown to 

/ \-l/2 / \-l/2 

be r(t) = ['2x\\t) at large t and with the solution of Eq.(|6|) at X'^ < 1 one obtains X{t) = Mbrf (2x||i] 
Upon substitution into Eq.(p7|) the running diffusion coefficient is obtained asymptotically as 



D{t-M,x\\,0.^)^{v\\)M^(2x\\t) (66) 



This is identical with the asymptotic behavior obtained from the quasilinear solution ([621). Thus, the stochastic 
generation of magnetic islands that appear at finite A_l does not affect either the asymptotic time-dependence of the 
running diffusion coefficient or its dependence on the parameters. 

There is however a significant effect of the nonlinear process of magnetic island generation but it appears to be 
localized in time. It can be found by determining the whole time evolution of the diffusion coefficient (|5^ ) using 
the computer code we have developed. The results are presented in Figures 1 and 2 compared to the solution (|6l|), 
( |6^ ) obtained for M <^ 1. One can see that at small and large times the diffusion coefficient is equal to Do{t). For 
intermediary times a transient decrease of D{t) appears. This is determined by the magnetic line trapping around 
stochastic magnetic islands, which is effective at times larger than the ffight time over the perpendicular correlation 
length Aj_, which in the unit considered here is Tfi = 1/M. As seen in Figures 1 and 2, the running diffusion 
coefficient has a maximum at Tfi and the Lagrangian velocity correlation becomes negative. Then the diffusion 
coefficient decreases due to the trapping of the magnetic lines which wind around stochastic island. This process 
is represented by the decorrelation trajectories corresponding to subensembles with large values of the parameter 
p ~ (jP /h which have performed many rotations around their paths (of small size) and their contribution cancels by 
mixing in the integrals in Eq.(|57|). Later in the evolution, another change of the sign of the Lagrangian correlation 
is observed at i = r^, the average return time for the parallel motion. At this moment D{t) has a minimum while 
DQ{t) has a maximum. It is determined by the parallel motion and more exactly by the collisions which force the 
particles to return on the magnetic lines. This is reflected in the decorrelation trajectories, which all evolve back on 
their paths in the perpendicular plane at t > t^. In the absence of the magnetic line trapping (quasilinear conditions) 
this leads to the decay of the running diffusion coefficient because the perpendicular displacement decreases in time 
and thus Do{t) decays aX t > Tr- The effect is inverse in the presence of stochastic magnetic islands. The backward 
motion produces first the un-mixing of the contribution of the trajectories that evolve on trapped magnetic lines. 
As time increases, the contributions of smaller and smaller magnetic islands are recovered in the Lagrangian velocity 
correlation. The effect of magnetic line trapping that produced the decay of D{t) in the interval (t/;, Tr) is washed out 
by the backward motion and D{t) recovers its value a.t t ^ Tfi. At this moment Tb, the correlation built-up time, D{t) 
has a maximum. A positive bump appears in the Lagrangian velocity correlation due to the trajectories unwinding 
around the magnetic islands. Finally, all decorrelation trajectories are "in phase" and approach the origin. This 
corresponds to the asymptotic regime in the evolution of the diffusion coefficient D{t) which is the same as for Do{t). 
Thus, the parallel coUisional motion eliminates asymptotically the nonlinearity determined by the x-dependence of 
the magnetic field ffuctuations. 

The above evolution of the diffusion appears whenever Tfi < Tr, and since Tfi — M^^ and ~ X\\ the condition 
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10"= 10° 10= J 10* 

FIG. 2. The running diffusion coefficient corresponding to tlie Lagrangian velocity correlations in Figure 1: Do{t) is the 
integral of Lo{t) and is given by Eq.(|62[) and D{t) is the integral of L(t) and shows the effect of the magnetic line trapping. 
The normalization constant is (A^J^) M. 

is Km > 1 which corresponds to magnetic hne trapping. When Tfi > Tr (or K„i ^ 1), the diffusion coefffcient is given 
byEq.dl). ^ 

We show in the next sections that this rather nontrivial evolution of the running diffusion coefficient leads to 
anomalous diffusion regimes when a decorrelation mechanism is present. 



V. DIFFUSIVE TRANSPORT INDUCED BY COLLISIONAL DECORRELATION 

We analyze in this section the effect of the cross-field collisional diffusion {x± 0) starting from the general 
solution (|56|)-(|5^). The stochastic collisional velocity ri±{t) in Eq.(^) moves the particles out of the magnetic lines 
and consequently it has a decorrelation effect leading to diffusive transport. This collisional motion determines 
a characteristic time, the perpendicular decorrelation time t±. It is defined by the condition that the collisional 
diffusion covers the perpendicular correlation length, 2xj_t± = A^, and in the units chosen here it is t± — (2x_l)~^. 
The stochastic magnetic field is considered here to be static (tc = oo) for a better understanding of the collisional 
decorrelation. 

As in the previous section, a stochastic magnetic field with small Kubo number Km that does not generate stochastic 
magnetic islands is first considered. We show analytically that the already known results are reproduced by the 
decorrelation trajectory method. Then the nonlinear case is analyzed and new anomalous diffusion regimes are found. 
They are determined by the non-linear interaction of the magnetic line trapping with the cross-field collisional diffusion. 

In 1979 Kadomtsev and Pogutse Q derived semi-qualitatively an approximation for the cross-field diffusion co- 
efficient. This approximation is essentially a weak-nonlinearity regime, in which the magnetic field fluctuations are 
non-chaotic. It will be shown that this diffusion coefficient is obtained from the general equations (|5^)-(|5|) provided 
that Tr < T± < Tfi. This condition is compatible with the relations found in |TT| where a detailed study of the 
diffusion regimes in stochastic magnetic fields for fusion plasmas is presented. In this conditions the X'^-dependence 
of the average velocity in Eqs.(j5^), (|5^) can be neglected and the equations for the decorrelation trajectories are 
( ^ ) corrected by a factor n'^{t) that multiplies the right hand side terms. This leads to the following form of the 
Lagrangian velocity correlation 

LKPit)^nl{t)Loit), (67) 

where Lo{t) is the subdiffusive Lagrangian velocity autocorrelation defined in Eq. (^l]). Because of the factor n^{t), 
the integral of Lxpit) no longer vanishes, and yields a finite diffusion coefficient, Dkp- It can be estimated analytically 
by using a step approximation of the function n±_ [t) 

nAt)^{ll Y>,\ (68) 
It then follows that the diffusion coefficient is approximated as: 

Dkp= dtLo{t) = - dtLoit) (69) 

Jo Jt± 
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FIG. 3. The asymptotic diffusion coefficient as a function of x±- Tlie total diffusion coefficient D (continuous lines) is 
compared witli the direct collisional contribution Xi_ (dotted line) and with the interaction term Dint (dashed lines) for two 
values of The normalization constant is X^i^, M = 10 and Tc = oo. 



because the integral of L(){t) from t — to infinity is zero. Using the very simple asymptotic form of Lalt) [obtained 
from Eq. ( |6l| ) for r > r^] , the integral can be calculated analytically and one obtains (going to dimensional quantities) 

Dkp = T^VWXI (70) 

which is the well-known Kadomtsev-Pogutse formula. 

When the time of flight Tfi is smaller than the decorrelation time t_i_, the space dependence of the magnetic 
field fluctuations cannot be neglected. It leads to stochastic magnetic islands. In the presence of a perpendicular 
collisional diffusivity the decorrelation trajectories obtained from Eq.(^) are not more closed curves. However, 
trajectory winding can still be observed for some range of the parameters that define the subensembles. This means 
that the process of generation of stochastic magnetic island and of magnetic line trapping still exists. Compared to 
the decorrelation trajectories obtained with x± — 0, these trajectories saturate faster and perform a smaller number 
of rotations. They still turn back at the maximum of the function T(i) which shows that the parallel trapping 
determined by the parallel collisional motion still exists. But due to the cross field collisional diffusion, these two 
trapping processes are only approximate or temporary. The perpendicular diffusion xa- produces a releasing effect both 
for perpendicular and parallel components of particle motion. The asymptotic values of the decorrelation trajectories 
are not concentrated in the origin (as for x± ~ 0) but spread in the {X, Y) plane. Consequently, a finite value of the 
asymptotic diffusion coefficient yields from Eq.(|57|). 



The asymptotic diffusion coefficient is determined from Eq.(57) using the numerical code we have developed. Some 
results are presented in Figure 3 where the asymptotic value of D(t) is represented as a function of x±- The two 
components Dint and x± are also represented. One can see that at small collisional diffusion Xi_ ^ Ij the non-linear 
interaction term largely dominates the collisional term while at large collisional diffusion x± ^ Ij the nonlinear term 
is only a correction to x±- Thus, the subdiffusive transport appearing at x± = is transformed by a small collisional 
cross field diffusion into a diffusive transport with a diffusion coefficient that can be several orders of magnitude 
larger than x±- The dependence of the diffusion coefficient on Xi_ is rather nontrivial. There is at very small Xj_ an 
increase of Z) up to a maximum which corresponds to t± = ti,. Then, at larger x±, the nonlinear interaction of the 
parallel and perpendicular trapping with the collisional decorrelation generates a strange transport regime, in which 
the effective diffusion coefficient decreases as the collisional diffusion x± increases. A minimum of D is obtained when 
X±_ determines a decorrelation time of the order of the return time of the parallel motion, t±_ = t^. At larger (when 
T_L < Tr), the nonlinear contribution Dint increases again with the increase of but this contribution begins to be 
comparable and eventually negligible compared to the collisional diffusion coefRcient. 

We note that the above results obtained with the decorrelation trajectory method are not similar with the heuristic 
estimation of the asymptotic diffusion coefficient of Rechester and Rosenbluth |^ . This is possibly due to the fact that 
the trapping of the magnetic lines, which is profoundly implied in the above results, is neglected in the estimation 0| 



and also in the more detailed calculations presented in |10|. This estimation is based on the process of exponential 



increase of the average distance between two magnetic lines in a chaotic magnetic field, represented by the Kolmogorov 
length. The estimation of this length taking into account the trapping of the magnetic lines should be necessary in 
order to compare the results. 
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FIG. 4. The asymptotic diffusion coefficient as a function of Tc for X\\ =0.1 (circles) and X|| = 1 (stars). The continuous 
lines represent the running diffusion coefficient as a function of t for the subdiffusive transport corresponding to static magnetic 
fields (tc — oo). The normalization constant is (A5_!^) M^, M = 10, Xi_ ~ 0. 



VI. DIFFUSIVE TRANSPORT IN TIME-DEPENDENT STOCHASTIC MAGNETIC FIELDS 



In a time-dependent stochastic magnetic field with finite Tc the configm'ation of the stochastic field b(x, z, t) changes, 
the magnetic lines move and consequently the perpendicular velocity of the particles is decorrelated leading to diffusive 
transport. We determine here the diffusion coefficient in such time-dependent fields in the limit of zero Larmor radius, 
stating from the general solution ( ^6|)-(p7|) . The effect of time variation of the stochastic magnetic field on the effective 
diffusion was previously studied in |]17|- but only for weak magnetic turbulence (if„ <C 1). We determine the 
effect of stochastic island generation appearing in stochastic magnetic fields at Km > 1. 

The decorrelation trajectories obtained from Eqs.(p3|), ( |5^ ) are in this case (finite Tc, — 0) situated on closed 
paths (except that for p = 0). A typical trajectory rotates on the corresponding path, then it stops and turns back. 
Its velocity decays progressively and eventually the trajectory stops somewhere on its path. This is the modification 
determined by the time variation of the magnetic field: all decorrelation trajectories stop at a time of the order Tc- 
Consequently, the running diffusion coefficient saturates. Depending on the relation between the decorrelation time 
Tc and the three characteristic times of this motion, r/;, r^, rt, (see Fig. 2) several diffusion regimes are obtained. In 
time-dependent magnetic fields, at i < Tj,, the time evolution of the diffusion coefficient is approximately the same 
with that obtained for Tc — > oo, and later, at t > Tc, D(t) saturates. Thus, the asymptotic diffusion coefficient can be 
evaluated as 



hm D{t;M,x\\,0,Tc)^D{Tc;M,xii,0,^). (71) 

using the running diffusion coefficient obtained in the static case. Thus it can be approximating by the value of the 
running diffusion coefficient for the subdiffusive case at t = Tc- Some results are presented in Figure 4 where the 
asymptotic diffusion coefficient obtained from Eq. ( |5^ ) for finite Tc is compared to the subdiffusive running diffusion 
coefficient represented in Figure 2. One can see that the approximation ( |7l| ) is rather good for all values of Tc- 

The following diffusion regimes can be observed in Figure 4, in the nonlinear conditions when stochastic magnetic 
islands are generated {tji < r^, or Km > 1). The quasilinear regime at small correlation times (tc < t//) with 
Do « M'^Tc is characterized by a fast time-variation which prevents trajectory trapping. At larger correlation 
times {Tfi < Tc < Tr) the stochastic magnetic islands can be generated before the stochastic magnetic field changes 
and the parallel motion is ballistic. In these conditions the diffusion regime is similar to that described in [|6) for 
the electrostatic turbulence: the diffusion coefficient decreases with the increase of Tc- A minimum of the diffusion 
coefficient appears at Tc = , followed by an anomalous increase determined by the interaction of the parallel trapping 
with the magnetic line trapping which generates correlation of the Lagrangian velocities. At very large correlation 
times (tc > n) the diffusion coefficient decreases as D « KmTc ^^^X\/'^ ■ We note that the regimes obtained for Tc < tji 
and for Tc > Tb are similar with those reported in [ pTf , |18|. But instead of the plateau found there at intermediary 
Tc, we obtain here a more complicated behavior. This is the effect of stochastic magnetic island generation: it leads 
to the decrease of the effective diffusion coefficient with the increase of Tc when the parallel motion is ballistic and, 
on the contrary, to the increase of D with the increase of Tc when the parallel motion is diffusive. 
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VII. CONCLUSIONS 



Wc have studied here the transport of coUisional particles in stochastic magnetic fields using the decorrelation tra- 
jectory method. We have derived analytical expressions for the running diffusion coefficient and for the Lagrangian 
velocity correlation in terms of a set of deterministic trajectories. They are defined in subensembles of the realizations 
of the stochastic field as solution of differential (Hamiltonian) equations that depend on the given Eulerian correlation 
of the stochastic potential. They are approximations of the subensemble average trajectories and represent the dy- 
namics of the decorrelation of the Lagrangian velocity. Since in general the equations for the decorrelation trajectories 
cannot be solved analytically, a computer code was developed for determining the running diffusion coefficient for 
arbitrary values of the four parameters of this problem and for given Eulerian correlation of the potential. 

We have shown that this rather complicated triple stochastic process is characterized by two kinds of trajectory 
trappings and contains two decorrelation mechanisms. The latter are produced by the coUisional cross field diflFusion 
X± and by the time variation of the stochastic magnetic field. 

One of the trapping processes concerns the parallel motion and is determined by collisions which constrain the 
particles to return in the already visited places with probability one. This parallel trapping leads to a subdiffusive 
transport in the absence of a decorrelation mechanism. This already known process is recovered by our method. 
The second kind of trapping concerns the magnetic lines which at Km > 1 wind around the extrema of the vector 
potential generating sclf-consistently magnetic islands. The effects of the magnetic line trapping in the presence of 
particle collisions is studied for the first time. Wc show that in the absence of a decorrelation mechanism, the stochastic 
magnetic islands determine a transitory decay of the running diffusion coefficient D{t) appearing at t in the interval 
{Tfi, Tr), i.e. before the parallel trapping is effective. The simultaneous action of both trapping processes determine 
a nonlinear built up of Lagrangian velocity correlation and eventually the parallel motion washes out the effect of the 
magnetic line trapping. Consequently, the asymptotic behavior of the running diffusion coefficient is exactly the same 
as in the quasilinear conditions when the stochastic magnetic field does not generate magnetic islands. 

The effect of the two decorrelation mechanisms is afterwards studied. We show that the effective diffusion coefficient 
and its dependence on the parameters results from a competition between the trapping and the decorrelation processes 
and more precisely from the temporal ordering of the characteristic times of these processes. Each one of the two 
decorrelation mechanisms leads to the already known diffusion laws when the stochastic magnetic islands arc not 
present (Km ^ !)• Their presence (at Km > 1) produces a complicated nonlinear interaction between the three 
stochastic processes which determines new scaling laws of the diffusion coefficient. They appear when the decorrelation 
time is longer than the flight time r/; but smaller than the correlation built up time Tt- The first condition ensures the 
magnetic islands generation and the second prevents the elimination of their trapping effect by the parallel coUisional 
motion. A particularly interesting regime is obtained for coUisional decorrelation and consists of an effective diffusion 
coefficient that decreases when the coUisional perpendicular diffusion increases (Fig. 3). 

This rather complex dependence of the diffusion coefficients on the plasma parameters can be used in experiments 
for controlling the transport. Even without changing the characteristics of the stochastic magnetic field, the diffusion 
coefficient can be strongly influenced by the parameters which describe particle collisions. A minimum of the diffusion 
coefficient was obtained for decorrelation times of the order of the average return time for the parallel motion. 
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